Formation of long-range spin distortions by a bound magnetic polaron 
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The structure of bound magnetic polarons in an antiferromagnetic matrix is studied in the frame- 
work of two-dimensional (2D) and three-dimensional (3D) Kondo-lattice models in the double ex- 
change limit {Jh ^ t). The conduction electron is bound by a nonmagnetic donor impurity and 
forms a ferromagnetic core of the size about the electron localization length (bound magnetic po- 
laron). We find that the magnetic polaron produces rather long-range extended spin distortions of 
the antiferromagnetic background around the core. In a wide range of distances, these distortions 
decay as l/r^ and 1/r^ in 2D and 3D cases, respectively. In addition, the magnetization of the 
core is smaller than its saturation value. Such a magnetic polaron state is favorable in energy in 
comparison to usually considered one (saturated core without extended distortions). 
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I. INTRODUCTION 

The physics of electronic phase separation is very pop- 
ular nowadays in connection with the strongly correlated 
electron systems, especially, in materials with the colossal 
magnetoresistance such as manganites^'^. The phase sep- 
aration manifests itself in other magnetic materials such 
as cobaltites^:, nickelates^ and also in low-dimensional 
magnets^. The nanoscale inhomogeneities were also 
widely discussed for heavy fermion compounds^, string 
oscillators in magnetic semiconductors^, electron-hole 
droplets in semiconducting nanostructures^, and high- 
Tc superconductoraSiifliii, in particular, in relation to 
the phase separation in the t — J modeli^ and to the 
formation of paramagnetic spin bags^- . 

The formation of small ferromagnetic (FM) metallic 
droplets (magnetic polarons or ferrous) in antiferromag- 
netic (AFM) insulating matrix was discussed beginning 
from works Refs. IT^flsl In application to manganites, 
this problem was addressed in a number of papers (see 
for the review Ref.P). Usually in manganites, the 'rigid' 
(well-defined) magnetic polarons with very rapidly de- 
creasing tails of magnetic distortion are considered. In 
other words, the intermediate region where the canting 
angle changes from (FM domain) to tt (AFM do- 
main) is narrow (of the order of interatomic distance d), 
so the radius of the magnetic polaron a is a well defined 
quantity. Moreover, the tails of magnetic distortions ex- 
ponentially decrease for the distances r > a (outside the 
ferron) . 

Long ago, de Gennes discussed the possibility of a 
slow decay of the AFM order distortions due to the 
magneto-dipole interaction^®. An attempt to get such 
a type of 'coated' magnetic polarons was made for a 
one-dimensional (ID) AFM chain in Refs-imi These 
calculations show that the characteristic length of the 



distorted spin surrounding can be much larger than the 
size of the trapping region. However, in ID, the 'coated' 
ferrons turn out to be metastable objects while 'rigid' 
magnetic polarons correspond to the ground state^^. A 
continuous 3D model of 'coated' magnetic polarons with 
the account taken of the tails of electron wave function 
in the AFM matrix (which appeared when we consider 
quantum nature of canting) was studied in Ref. l2ri Un- 
fortunately, even in this more advanced approach, the 
tails of the magnetic distortions exponentially decrease 
outside the ferron. However, in Ref.|20, the gradient con- 
tribution to the magnetic energy proportional to (Vz^)^ 
was not included to the variational procedure to reduce 
the order of the differential equations under study. In the 
present paper, we generalize the ID model considered in 
Ref. to the 2D and 3D cases. In contrast to the ap- 
proach of Ref. I23, we take into account the magnetic 
anisotropy, discrete nature of the lattice, and the gra- 
dient terms. Appropriate account of the gradient terms 
leads to a slower decrease of magnetic distortions outside 
the ferron. 

In our model, we consider an electron, which is bound 
at a donor impurity by the Coulomb attractive poten- 
tial. The Coulomb potential V is assumed to be strong 
in comparison to the other relevant interactions. Namely 
F ~ J_ff ^ i ^ Jdd, K, where Jh is the Hund's rule cou- 
pling, Jdd is the AFM exchange interaction, t is a hopping 
integral, if is a constant of the magnetic anisotropy. In 
this range of parameters, the radius of the electron lo- 
calization is of the order of interatomic distance d. The 
magnetic distortions outside the electron localization re- 
gion (ferron core) are shown to decrease with the distance 
r as l/r"' and in the 3D and 2D cases, respectively. 
In contrast to the ID model^^i, these 'coated' ferrons cor- 
respond to the ground state of the system. 

Note that the 2D case widely discussed in our pa- 
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per could be rather relevant to experiment since a lot 
of magnetic oxides including manganites and cobaltites 
have layered magnetic structure and a 2D structure is 
simply a limiting case of the layered system. Actually, 
as it was shown by neutron diffraction studiesSLSSiSiSi, 
the ferrous in low-doped manganites could be 'platelets' 
characteristic for 2D rather than spherical droplets char- 
acteristic for 3D. 

At the end of the paper, we find the applicability range 
of the proposed approach by taking into account the tails 
of the electron wave function outside the ferron. Our 
results enrich the concept of microscopic phase separation 
for the localized magnetic polarons. 



II. THE MODEL HAMILTONIAN 

We consider 2D or 3D cubic lattice of antiferromag- 
netically coupled local spins Sn treated as classical vec- 
tors. We assume that the crystal has uniaxial magnetic 
anisotropy, with x axis being the easy axis. The nonmag- 
netic donor impurities are placed in the center of some 
unit cells of the lattice. It is assumed that the doping con- 
centration is small enough and therefore we can consider 
an isolated impurity and restrict ourselves to a single- 
electron problem. The Hamiltonian of such a system has 
the following form: 



H — Hcl + Jdd ^ [SnSm + S"^] 



(1) 



(nm) 



-K 



(niTi)CT 



|n- no| 



(2) 



where aL 



are the creation and annihilation opera- 



site n, a are Pauli matrices, and the symbol (...) denotes 
the summation over nearest neighbors. The second and 
third terms in Eq. are the AFM exchange between 
local spins and the magnetic anisotropy energy, respec- 
tively. The first two terms in Hc\ describe the kinetic 
energy of conduction electrons, and the Hund's rule cou- 
pling between the conduction electrons and the localized 
spins. The last term in Eq. Q describes the Coulomb in- 
teraction between electrons and the impurity ion placed 
in the middle of the main diagonal of the unit cell (e.g., 
in the 3D case, between lattice cites n = {0, 0, 0} and 
n = {1, 1, 1}, that is, no = {1/2, 1/2, 1/2}). 

Below we consider the range of parameters Jh ^ t ^ 
Jdd ^ K' characteristic of manganites. In the limit 
Jh ^ oo, the spin of electron at site n is parallel to 
the local spin Sn = S'{sin6'n cost^n, sin 0^ sin (/3n , cos^n}. 
This implies the transformation of Una to operators 
with spin projection -1-1/2 onto the direction of Sn (spin- 
less fermions) and also to the transformation of hopping 
amplitudes: t — + t cos(i^nm/2)e^*'^""', where j/nm is the 
angle between Sn and Sni 

cos l/nm = COS 9n COS 6*^ + Siu ^n siu 9^ C0s((^n — <^m), 

and 



arg 



COS ■ 



■ cos ■ 



-I- sin — sin ■ 



tors for the conduction electron with spin projection a at 



is the Berry phase2S*2&. 

At low doping, the electron is bound by the impurity 
electrostatic potential. In the limit of strong electron- 
impurity coupling V ^ CO, the electron wave function 
will be nonzero only at sites nearest to the impurity. Let 
us focus first on the 2D case. The three-dimensional case 
is considered in the similar way. Supposing that ^'n 7^ 
only for ni = {1, 1}, 112 = {1, 0}, ns = {0, 0}, and 
n4 = {0, 1}, and following the standard diagonalization 
procedure, we find for smallest eigenvalue of the electron 
energy 



+ V [(Cl2 - C34)2 + (C23 + CAlY] [(Ci2 -|- €34)^ -|- (C23 - 041)^] - I6C12C23C34C41 siU 



1/2 



(3) 



where the constant term — JhS/2 — V\/2 is omitted. In 
expression (O, = cos(i^ninj /2), tOij = ujmnj, and 
Aw = W12 -I- 1^23 + ^34 -I- a;4i. Note that in the general case 
we have Alu ^ 0. The electron energy E^i has a minimum 



when all spins Sn; are parallel to each other. Thus, we 
have a bound magnetic polaron state: a ferromagnetic 
core embedded in the antiferromagnetic matrix. 
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FIG. 1: 'Coated' magnetic polaron before (left figure) and 
after (right figure) transformation of angles in one sublattice. 
The magnetic structure is calculated by solving Eqs. ® at 
t/J = 50, X — 5 X 10"'^. At this values of parameters, (po = 
85°. 



III. MAGNETIC STRUCTURE IN THE 2D 
CASE 



In this section, we calculate in detail the magnetic 
structure of bound magnetic polarons for the 2D lattice. 
The problem can be solved by minimization of the total 
energy with respect to angles 9n and (p^. It is conve- 
nient to make first the following transformation. In the 
absence of magnetic polaron, the lattice consists of two 
magnetic sublattices with antiparallel magnetizations. In 
one sublattice, we perform the transformation of the an- 
gles 9n ^ T^ — On and ipn mod(7r-|-(/Jn, Stt). As a result, 
an AFM order becomes FM, and vice versa. Such a trans- 
formation allows us to work with continuously changing 
orientation of spins outside the ferron core. The total 
energy then reads 

E = Ecl + J ^ [1- cos J/nm] 
(nm) 

-K [sin^ On cos^ (p„ - l] , (4) 



where J = JddS'^, and K = K'S'^. Note that after the 
transformation, we should replace in Eq. Q — > Sij = 
sin(i/ninj /2) and Aw — > —Auj. 

Let the lattice plane be xy plane. In the range of pa- 
rameters under study, the trivial solution of the problem 
Eq. I0J corresponds to the case On ~ 7i'/2, = for 
n ^ rii, and = i^na = 0, (pn2 = y^n^ = tt. So, we 
have a bound magnetic polaron with completely polar- 
ized spins embedded in purely antiferromagnetic back- 
ground. The total magnetic moment of such a polaron is 
parallel to the easy axis. We refer this trivial solution as 
a 'bare' magnetic polaron. Its energy is = —2t+ 16 J. 

There is another solution corresponding to the mag- 
netic polaron state with magnetic moment perpendicu- 
lar to the easy axis: a 'coated' magnetic polaron. Let 
us assume again that all spins in the lattice lie in the xy 
plane, that is, = 7r/2. It is clear from symmetry that 
angles ipi = ip^- for the local spins inside the magnetic 



polaron should satisfy the conditions 



(5) 



where Q<ifQ<TT/2 (see Fig. Minimizing the total 
energy Eq. I^J with respect to ip-n = ^n^,ny, taking into 
account relation ^ for Ec\ and the fact that Auj ~ for 
a planar configuration of spins, we find the following set 
of nonlinear equations 



sin((^n+A - Vn) - ^ sin 2(p„ 



(6) 



where k — 2K/ J, Snm is the Kronecker symbol, and A 
takes values {±1, 0}, {0, ±1}. 

Equations ((HJ with conditions ^ are solved numeri- 
cally for the cluster containing 40 x 40 sites. The fur- 
ther growth of the number of sites in the cluster does 
not change the obtained results. The initial angle ipo is 
also found. The calculated magnetic structure is shown 
in Fig. n We see from this figure that, in contrast to 
the 'bare' ferron, the 'coated' magnetic polaron produces 
spin distortions of AFM background outside the region of 
electron localization. In addition, the 'coated' magnetic 
polaron has a magnetization lower than its saturation 
value {ipo < 7r/2). Moreover, the 'coat' has a magnetic 
moment opposite to that of the core (sites 1-4 in Fig.^). 

In order to get analytical estimations for the spa- 
tial distribution of the spin distortions, we find an ap- 
proximate solution to Eqs. © in the continuum limit. 
Namely, angles (/?„ are treated as values of continuous 
function (p{r) at points r = n — Hq. Assuming that out- 
side the magnetic polaron the following condition is met 
|(/5n_i_A — Vnl ^ 1, we can expand ip{r + A) in the Taylor 
series up to the second order in A 

ip{v + A) « ^(r) + A'^dMr) + ^A" A^9,9^^(r) . (7) 

Substituting this expansion into Eq. ®, we find that 
function ip{r) outside the magnetic polaron should satisfy 
the 2D sine-Gordon equation 



Aip — — sin 2ip — . 



(8) 



In the range of parameters under study, K <^ J, that is 
>c 1, we can linearize this equation. As a result, we 
obtain 



Aip — xip = . 



(9) 



This equation should be solved with the boundary con- 
ditions at infinity (p{'r)\r^oc-^ and with some bound- 
ary conditions at the surface of the magnetic polaron. 
We model the magnetic polaron by a circle of radius 
a — 1/ \f2 (in the units of lattice constant d) and choose 
the Dirichlet boundary conditions 



(10) 
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FIG. 2: Angles (/Jn.n (circles) calculated by solving numeri- 
cally Eqs. ® at tjj = 50, >ir = 5 X 10~^ (ro ~ 4.5). The curve 
corresponds to the function <^(n — 1/2, n — 1/2), see Eq. II14|I . 
The initial angle ipo in Eq. 1141 is used as a fitting parameter. 
The range of distances where the quantum fluctuations can 
play an important role at S = 2 is shown by the dashed line. 



FIG. 3: ip{x,y) vs ^ = arctan(x/?/) at x + y — N for 
= 3, 4, 5, 7 from the top to bottom calculated accord- 
ing to Eq. I|14|l (dashed curves). Points correspond to the 
numerical calculations of angles Lp„^M-n+i (Cn = arctan[(n — 
l/2)/(A - n + 1/2)]). Parameters of the model are t/J = 50, 
M^5x 10"^ 



where we introduce polar coordinates (r, C) in the xy 
plane. The function (^(C) can be found in the follow- 
ing way. Note, that <^(C) should satisfy the symmetry 
conditions lO at points Ci = ""(Si — l)/4 

t = l..A. (11) 

Since the function (p(C) is a periodic one, it can be ex- 
panded in the Fourier series 

oo 

(p{C) = ^ [a,„ COS mC + 6m sin mC] . (12) 

m=0 

In Eq. (|12|l we neglect the terms with m > 2 which allows 
us to keep the minimum number of terms to satisfy the 
conditions (|ll|l . It follows from Eq. Hll|l that qq = ai = 
02 = 6i = 0. Finally, we obtain 

(p(C) =(posin2C. (13) 

The solution to Eq. © with boundary condition ^W^i 

is 

= irrTT^2 (-) sin 2C, ro = ^ , (14) 
K2(a/ro) \roJ 

where K2{x) is the Macdonald function. Within the 
range |r| < tq, (p{r) behaves as a'^/r'^, whereas at large 
distances, it decreases exponentially (fir) cx exp(— r/rp). 
The function 95 (r) at ^ = 7r/4 and the numerical results 
for ifn^n are plotted in Fig. |21 In order to compare the 
angular dependence of the approximate solution with the 
numerical results, we plot in Fig.|31the values oiip{x, y) at 
lines X + y = N depending on the angle ^ — arctan(x/?/) 
at different N. The values of ifn,N-n+i are also shown 
in this figure. We see from these figures that, except 
for the small region near the magnetic polaron, the con- 
tinuous function Eq. (|14|) is a good interpolation of the 
function (/5„^.„^ at the discrete lattice. Note that at small 
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FIG. 4: The energy difference AE between 'coated' and 'bare' 
ferrous ys — 2K/J calculated a,t t/J = 50. Solid curve cor- 
responds to the numerical calculations, whereas dashed curve 
is calculated using Eq. 1161 . 



r, the continuum approximation fails because of the large 
derivatives of ip(r). In addition, we cannot neglect the 
nonlinearity of differential equation (jS)) in this region. 

Let us calculate now the energy of 'coated' magnetic 
polaron, i?p, in comparison to the 'bare' one = 
—2t+ 16 J. The energy of the system is given by Eq. 
where 9^ = 7r/2 and is a solution to Eq. 0. In 
the continuum approximation, we expand cos Vnm and 
cos^ (Pn in the Taylor series up to the second order outside 
the ferron core. Changing the summation by integration 
over r, we find for the last two terms in Eq. 

8 J (1 -I- 1) sin2 ^o + ^J d\ [{V^f + V] • (15) 

r|>a 

In Eq. H15|l , the first term comes from the summation over 
spins in the core of the magnetic polaron. Substituting 
solution (|14|l into Eq. H15|l and performing the Integra- 
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tion, we find 



where 



^) sinVo + ^/(aV^)-2isin(^o, (16) 



Hpo) = 27r 



1 



(17) 



The initial angle ipo is found by minimization of en- 
ergy ()16fl . This value turns out to be slightly differ- 
ent from that found numerically. The dependence of 
AE = Ep — Ep on X is shown in Fig.0J We see from this 
figure that AE < at not too high values of anisotropy 
constant K. So, the magnetic polaron with extended 
spin distortions can be more favorable than the 'bare' 
one. Note also that formula H16|l overestimates the en- 
ergy of the 'coated' magnetic polaron. 



IV. MAGNETIC POLARON IN 3D LATTICE 

In this section, we calculate the magnetic structure 
of the 'coated' magnetic polaron in the 3D case using 
the approach described above. In the V oo limit, a 
magnetic polaron includes 8 sites, that is, ^ only 
for m = {1, 1, 1}, n2 = {1, 0, 1}, na = {0, 0, 1}, n4 - 
{0, 1, 1}, ns = {0, 1, 0}, ne = {1, 1, 0}, = {1, 0, 0}, 
ng — {0, 0, 0}. The electron energy is found then by the 
diagonalization of the corresponding 8x8 matrix. As in 
previous section, we assume that all spins in the lattice lie 
in the xy plane. The symmetry conditions then become 



(18) 



Minimizing total energy Q with respect to (^n , we get 
a system of equations similar to Eqs. © with the only 
difference that the factor at the sum in the right-hand 
side is equal now to 3t/8J and A takes values {±1, 0, 0}, 
{0, ±1, 0}, and {0, 0, ±1}. We solve this system of equa- 
tions numerically for the cluster with 30 x 30 x 30 sites. 

In the continuum approximation, the magnetic struc- 
ture is found in the same manner as in the 2D case. 
We consider a spherical magnetic polaron of radius a = 
\/3/2. The function (p{r) satisfies Eq. with the bound- 
ary conditions at |r| — a 



f{r)\r=a = ^(t?,C) , 



(19) 



where r, and C are spherical coordinates and the func- 
tion (^(??, C) should satisfy the symmetry conditions l|18|) 
at angles Si, Ci corresponding to the directions of n.; — Uq 
{i = 1...8). Expanding (p{'d,C) in series of spherical 
fimctions Yimi'&iC) ^-nd retaining only the terms with 
the smallest 1, one finds 



mo 



36n 



¥'0V^[Y3,+2(^,C)-Y3.-2(^,C)] 

'27 2 
(fio \l — sin d cos d sin . 



(20) 
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FIG. 5: The dependence of <^(r) on the distance r in the 
direction {1, 1, 1} calculated according to Eq. 12111 (solid and 
dashed curve). Circles represent the numerical results for 
at t/J — 100, X = 5 X 10 ^. The initial angle ipo 
in Eq. 1211 is used as a fitting parameter. The range of dis- 
tances where the quantum fluctuations can play an important 
role at S = 2 is shown by the dashed line. In the inset, the 
dependence of AE on is shown. The energy difference AE 
is calculated numerically at t/J = 100. 



The solution to Eq. @ with the boundary conditions l|19l) 



(f{x,y,z) = ifo 



(21) 



where again tq = 1/ \/>c, and 



^ + ^)-r- (22) 
5 15 / p4 y ' 



It follows from Eqs. (|21|l and (|22|l that Lp{r) decreases 
at r < tq as a^/r^. At distances r exceeding character- 
istic radius tq, the spin distortions decay exponentially, 
</j(r) cx exp(— r/ro) as in the 2D case. The dependence 
of (^(r) on distance r in the rii = {1, 1, 1} direction is 
shown in Fig. El The numerical results for ipn,n,n are also 
presented in this figure. 

The energy of the magnetic polaron, Ep, is found in 
the same way as in the previous section. Note that in 
the 3D case, the energy of a 'bare' magnetic polaron is 



E" 
p 



-3t + 48 J. The dependence of AE = Ep 



E° 
p 



on is shown in the inset to Fig. |S| We see that in the 
3D case, the situation is similar to that considered in the 
previous section. At not too high anisotropy constant, 
the 'coated' magnetic polaron is favorable in energy in 
comparison to the 'bare' one. 

Note, however, that the two-sublattice AFM state con- 
sidered here is not eigenstate of the system, and quan- 
tum fluctuations play an important role. Indeed, the 
fluctuation magnitude is about SS/S — 0.078/5 for the 
isotropic 3D Heisenberg modelSi. This value becomes of 
the order of the angle (/) at r = 2 (second coordination 
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sphere). However, we consider the situation with uniax- 
ial anisotropy, which drastically reduces quantum fluc- 
tuations stabilizing the two-sublattice AFM state. For 
example, for manganites, the characteristic value of the 
anisotropy parameter K — 0.05 J (see Ref . l28|) . At such 
anisotropy, the fluctuation magnitude decreases by about 
30%. This means that we can use our results up to r = 3. 
At the concentrations of the impurities higher than 1% 
the distance between impurities is less than 6 lattice con- 
stants and the behavior of the distortions at r > 3 is 
not of real physical importance. The similar arguments 
are valid also for the 2D case. Here, we can find that 
5S/S = 0.092/5' at = 0.05J (note that the magnetic 
anisotropy stabilizes the long-range magnetic order in the 
2D Heisenberg model). In 2D, the distortions of AFM or- 
der decay slower than in the 3D case. Thus, the quantum 
fluctuations become comparable to the angle at r ~ 4. 
The range of distances where the quantum fluctuations 
can play an important role at 5 = 2 is shown by the 
dashed line in Figs. |21 and O 



V. CONCLUSIONS 

We studied the structure of bound magnetic polarons 
in 2D and 3D antiferromagnets doped by nonmagnetic 
donor impurities. The doping concentration was as- 
sumed to be small enough, and the system is far from 
the insulator-metal transition. It was also assumed that 
the conduction electron is bound by electrostatic poten- 
tial of impurity with the localization length a close to 
the lattice constant d. We found two possible types of 
magnetic polarons, 'bare' and 'coated'. At not too high 
constant of magnetic anisotropy K J, the 'coated' 
polarons correspond to a lower energy. Such a 'coated' 
ferron consists of a ferromagnetic core of the size a and 
extended spin distortions of the AFM matrix aro und it . 
The characteristic radius of distortions is tq = d^fjflK. 
Within the a < r < tq range, these spin distortions be- 
have as a^/r^ in the 2D case and as jr^ in the 3D 
case. At larger distances, they decay exponentially as 
exp(— r/ro). Note that in real systems, the constant of 
magnetic anisotropy K is small enough in comparison to 
exchange integral J and the characteristic radius ro can 
be rather large. 

The calculations have been performed in the case of 
zero temperature T = 0. At finite temperatures, the reg- 
ular picture of spin distortions could be, in principle, de- 
stroyed, owing to the spin-wave excitations. However, in 
the case of uniaxial anisotropy, the spin-wave spectrum 
has a gap of the order of \J2KJ. So, the characteris- 
tic temperature, below which our results (in particular, 
those presented in Figs. [5] and [3) are applicable, is of the 
same order. For example, in manganites, a typical value 
of this temperature is 5 K. 

Note that in the ID case, the 'coated' ferron was found 
to be a metastable objeclii^ in contrast to the 2D and 3D 
cases considered here. The cause of such a difference is 



the following. In the ID case, the 'bare' magnetic polaron 
has zero surface energy since it does not disturb the AFM 
ordering at its border. For I? > 1, the surface energy is 
nonzero due to the geometry of the problem. For the 
'coated' magnetic polaron, the surface energy is nonzero 
at any dimension. However, the extended 'coat' reduces 
its value. 

Our results were obtained in the limit of strong 
electron-impurity coupling, y — > cxd. In this case, the 
wave function \l/n of the conduction electron is nonzero 
only at the sites nearest to the impurity, and one can 
calculate the electron energy i^ei exactly. At finite F, 
the wave function extends over larger distances and one 
should calculate the magnetic structure simultaneously 
with ^fn. Note, that even at = 0, the magnetic po- 
laron is a stable object in a wide range of parameters of 
the model. This is the so called self-trapped magnetic 
polaron, which exists due to trapping of a charge carrier 
in the potential well of ferromagnetically oriented local 
spins The radius of such a polaron was estimated as 
a ~ d(tl 3)^1^ and at large tj 3 \\ can be of the order of 
several lattice constants^ 

Let us find the range of validity for the assumption 
V ^ oo used in this paper. For this purpose, we rewrite 
the electron energy i?oi in terms of the wave function. 
Performing the transformation of angles described in Sec- 
tion [ffll and using the substitution \l/„ = e*'''"/^V'n, we 
can write at 0„ 7r/2 

(nm) ^ ^ n I 

Minimizing the total energy with respect to angles ip-a^ 
we obtain the set of equations similar to Eqs. © with 
the right-hand side in the form 

t ^ In, (^„^;^^) . (23) 

In the limit of strong electron-impurity coupling, we can 
assume that the wave function is hydrogen-like, that is, 
IV^nl K, Np ' exp[— (d|n — no| — a)/a*g], where A^p is the 
number of sites inside the magnetic polaron and a|j is the 
effective Bohr radius. The latter one can be estimated 
as a*g = h^e/m*e'^, where e is the permittivity, and m* 
is the effective electron mass, which can be found from 
the relation /2m* = td^. Taking into account these 
estimations, we find that one can neglect the right-hand 
side in Eq. for n 7^ if 



In the 3D case, we have A^p — 8, and the approximation 
F CX3 is valid for V = e^/ed down to Vc = 2t\n{t/163). 
For realistic values of parameters V is usually of the or- 
der of Vc- So, it is of interest to perform a more detailed 
analysis of the problem at finite values of V similar to 
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that discussed in Ref. |^ for the ID case. Nevertheless, 
we hope that our results concerning the extended spin 
distortions are valid also a.t V < Vc, and, consequently, 
for magnetic polarons with larger sizes a. Indeed, in the 
continuum approximation, the dependence of the mag- 
netic structure on a is not so important. Moreover, the 
energy of the magnetic polaron Ep formally depends on 
a only through the term J(/7g/(a-\/x)/2 ~ J coming from 
the spatial integration (see Eq. HlGfl ). Since <C 1, the 
energy difference AE is negative up to a ~ rp ^ d. 

In our paper, we studied the isolated ferrous. However, 
as it was already mentioned above, the substantial over- 
lap of the extended distortions can arise even at very low 
doping (about 1%). As a result, the magnetic moments 



of the ferrons could be ordered. In particular, according 
to Ref. Is^ in such a situation, there can exist a long- 
range ferromagnetic order with the magnetic moment far 
from the saturation value. 
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